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Abstract 
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1 INTRODUCTION 



It is plausible that, in high-energy collisions, one occasionally produces a 
bubble of classical pion field, randomly oriented in isospace [|l], 0, §]. There 



is a growing literature on this subject 0-|[T^. All these admittedly very 
speculative papers are not purely esotheric, however. A transcient production 
of a " disoriented chiral condensate" could be an explanation for the heretofore 



mysterious Centauro events reported by cosmic ray people [jTT|. Indeed, it 
is easy to convince oneself that the fraction r of neutral pions produced in 
the decay of a bubble of disoriented chiral condensate should be distributed 
according to the probability law 

dP{r) = ^dr (1) 

Hence, the probability of events with small r is unusually high compared to 
naive statistical expectations. 

The specific scenarios proposed in the theoretical papers quoted above 
are more complementary than equivalent. Hopefully new data will select 
the right one. For the time being, it seems reasonable to develop all the 
a priori plausible scenarios, checking their self-consistency and looking for 
implications. Thus, the aim of this paper is to further develop the ideas 
put forward in ref. 0. We follow a strategy similar to that adopted there, 
sacrificing generality in favour of simplicity. 

The plan of the paper is the following: In sect. 2 we recall some of 
the results obtained in ref. |0| in 1-1-1 dimensions and we discuss a possible 
extension to 1-1-3 dimensions. We argue that the disoriented chiral condensate 
should be considered as an open quantum system. In sect. 3 we formulate 
the analogy between our problem and the kinetics of a laser under threshold. 
In sect. 4 we establish a hierarchy of time scales relevant to the evolution of 
the condensate. Sect. 5 contains the discussion. 



2 THE ORIGINAL MODEL AND GENE- 
RALISATIONS 

The starting idea is to use the classical equations of motion derived from an 
effective chiral Lagrangian to calculate the intensity of soft-pion radiation in a 
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high-energy heavy-ion colhsion. In ref. the lowest order chiral Lagrangian, 
viz. that of the non-hnear a model has been employed: 

^ = \fl + (5^)1 (2) 

and 

+ = 1 (3) 

The classical equations of motion are summarized in the following current 
conservation equations 



9V = 

(9A = (4) 

where and are the Noether isovector and isoaxial currents associated 
with the global SU{2) x SU{2) symmetry of the model. These equations 
must be supplemented by appropriate boundary conditions. Inspired by 
Heisenberg's old papers |T^ on semi-classical description of multi-particle 
production, we have adopted the following idealization: At time t = the 
whole energy of the collision is localized within an infinitesimally thin slab 
with infinite transverse extent (instead of a pancake shaped region). The 
symmetry of the problem then implies that the pion field depends only on 
the invariant s = t"^ — x"^, where x is the longitudinal coordinate. Thus the 
problem becomes effectively 1 + 1 dimensional and the currrents take the 
form 



A, = b d^(t> (5) 
where a, h are integration constants satisfying a ■ b = and 

9^0 = 0, s > (6) 

Hence, the problem is reduced to an Abelian one. The solution of @ corre- 
sponding to currents non- vanishing within the light-cone is (j){s) = log (s/sq). 
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Projecting the pion field on the axes of the triad a,b, c = a x b and after 
some algebra one finds for s > 



TTa = 

TTb = — sin (/€</)) 

TTc = (o/k) cos (7) 

with K = y/aFT¥. Of course, the classical solution breaks the symmetry 
of the theory: within the light-cone the pion field rotates in isospace in a 
plane perpendicular to a randomly chosen vector a. Outside the light-cone 
a = 1. The solution is sigular for s = . Hence, it describes radiation 
by sources living on the light-cone. For t = const, the pion field fills the 
space between these two sources and oscillates violently as one approaches 
them. One then enters the region where is stored almost all the energy of the 
collision and where the gradient expansion involved in the derivation of the 
effective Lagrangian from QCD stops making sense. Close to the sources the 
a model (|l]) cannot be trusted and, in our opinion, its use is not justified. 

Thus, the picture in 1-1-1 dimensions is essentially the same as in 0: there 
is an inner region, where resides a disoriented chiral condensate, separated 
by a potential barrier from the outer physical vacuum. All this comes out 
quite naturally. Unfortunately, the generalisation to 1 -|- 3 dimensions is not 
unique. One can assume as in refs. [|[ H that somehow the same topology 
is created also in this case. This occurs perhaps in events with exceptionally 
high multiplicity. However, in this paper we shall keep working with the 



boost-invariant scenario [0, |T3| as in 0], postulating that the sources are 
localized in those space regions where energy density is high in a generic 
heavy- ion collision (i.e. essentially within the two receding pancakes). We 
shall see soon that this assumption is less innocent than it might appear and 
leads to difficulties in the purely classical framework. The way out will be 
to consider the disoriented chiral condensate as an open quantum system, 
which it actually is. 

Before going further we need to consider for a while the modifications 
of the results obtained in 1 -|- 1 dimensions, when the initial system has a 
finite space extent e. As in we shall calculate the total radiated energy 
at some large time t, but we shall not include in the space integration the 
regions where reside the sources. Technically, in calculating the spectrum. 
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we Fourier transform the fields multiplied by a cutt-off function f{t, x) which 
equals 1 within the light-cone, except in the regions |x| < |t| — e where it 
vanishes rapidly. In analogy to we are led to calculate the integral of the 
type 

I = [ dxf{t - x)exp[ikx + a In (t^ - x^)], k > (8) 
Jo 

(for simplicity of writing we consider the integration over x > 0). Change 
the integration variable to z = t — x and consider the integration over the 
closed contour made up of three parts : < Re^; < t, — | < arg2; < and 
—t< Imz < 0. For large t one has 

/ ^ _,(_2^t)"e^'=* f rfi//(-^y)e- V[l + 0{y/t)] (9) 
Jo 

When / = 1 one can obtain from (j^) the leading term of the large-t expansion 
of the exact expression quoted in [Q (for t ^ oo the integral converges 
towards r(l + a)/k^~^°' ). When e > (and f{z) is smooth enough, e.g. 
f{z) ~ exp{—e/z)) the contribution to the integral in (^) from the region 
< ?/ < e is vanishing and there appears a factor exp{—e\k\) in the spectrum. 
This agrees with the physical intuition: only modes with wavelength larger 
than the size of the emitter are effectively produced. One could arrive at the 
same conclusion checking that the stationary point of the exponential in the 
integrand of (|) is at a distance of order 0(1/ A;) from the light-cone and is 
therefore outside of the cut-off region provided k is small enough. 

In 1 + 3 dimensions we have been unable to find the most general solu- 
tion of eqs. (^) for some physically reasonable set of boundary conditions. 
However, if one accepts to limit oneself to solutions of the form (|^), then 
the problem again becomes Abelian and the field equations reduce to (|^) 
(outside of sources). We have studied various solutions of the wave equation 
(^. The upshot of this study is almost evident. Therefore, we shall not enter 
into details, limiting ourselves to pointing out what is truly relevant. The 
solutions of interest are those with cylindrical symmetry : (p = (j){s, r), where 
r is the transverse coordinate. The initial transverse radius of the source is 
denoted by R. It is not obvious whether R should be assumed to be roughly 
equal to the nuclear radius of the colliding heavy ions or should rather be 
set to a smaller value. It depends on how much coherence one expects to 
have in the transverse direction and this in turn depends on the mechanism 
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triggering the creation of the disoriented chiral condensate. Most hkely, R 
fluctuates from event to event and there is a penalty for R being large (in 
nuclear units) . The generic behaviour of is 

0~logs, v^,r <i? (10) 

and 



~ s-\ y/s » R,r (11) 

Eq. (10) simply states that the one-dimensional idealization of ref. is 
roughly realistic as long as the sources are separated by a distance less than 



R. Eq. (|TT]) asserts that very far from a finite size source the field becomes 
that of a point emitter. Notice, that as times goes on the bulk of the field 
stays at a finite invariant distance from the sources. In the analogue of eq. 
(H) the stationary point is at a distance from the source, which inevitably 
becomes less than e provided t is large enough and whatever is k. In other 
words for t oo the field collapses on the sources and there is no radiation. 
A straightforward generalisation of the discussion of ref. 0] simply does not 
work. Of course, it is completely irrealistic to imagine that the two sources 
move indefinitely intact and with the velocity close to that of light. However, 
we do not see any reason why they would switch off suddently at t ~ i?, before 
the radiation enters into the three-dimensional regime. 

Trying to understand the physics of this problem, one should notice that 
the soft-pion radiation we are considering differs in (at least) one essential 
aspect from the familiar bremsstruhlung from a classical current discussed 
in any elementary textbook of quantum electrodynamics. The point is that 
the pion system is permanently in contact with the debris of the colliding 
nuclei. Thus the pion system is an open one and is subject to decoherence. 
We shall see in the next sections under which conditions the latter can occur 
before the one- dimensional expansion is over. 



3 THE DAMPED OSCILLATOR ANALOGY 

We shall focus in this section on the decoherence problem. Because the 
subject is perhaps not very familiar to all the potential readers of this paper, 
and also in order to spell out the numerous simplifications we are going to 
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adopt, we shall briefly sketch certain derivations. For more details and for 
references to the original works the reader should consult, for example, refs. 

PI, mm. 



The idea is to consider the disoriented chiral condensate as a " system" in 
contact with thermalized hadronic matter produced in a heavy-ions collision 
(the "bath"). We shall represent the "bath" by a collective variable B, which 
we shall couple linearly to the "system". One can write B = (B) + 6B . 
The classical component of B will yield sources responsible for the soft-pion 
radiation discussed in P]. The quantum fluctuations of 6B will tend to 
break the coherence of the radiation. Actually, in order to avoid technical 
complications we shall consider an Abelian toy model in 1+1 dimensions. The 
non-linear a model is reducible to this model under circumstances mentioned 
in the preceding section, but in the classical limit only. However, we hope that 
the toy model is sufficient to produce rough estimates of the characteristic 
time scales we are interested in. 

The starting point is the von Neumann-Liouville equation satisfied by the 
density matrix, denoted by W: 

W{t) = -i[H,W{t)] (12) 

where H is the Hamiltonian describing a free scalar field coupled linearly to 
the collective variable representing the "bath". We assume that the coupling 
(in the Lagrangian) is a gradient one, viz. ~ d^(f)d'^B, so that only a rapidly 
varying "system" field is strongly coupled. 

We partition the rapidity space into cells of extent 6y ^ 1. We denote 
by the operator creating a quantum of our scalar field within a given cell 
and we write 

H = ua^a + Hbath + V (13) 

with the interaction 

V = {a + a^)j{t) + X{uj)[ab'^ + a^b] (14) 

Here, u is the average energy within a cell while j{t) and b are the appropriate 
projections of the classical current j ~ (d'^B) and of the quantum "bath" 
field 6B, respectively. We assume that Hbath does not couple distinct cells. 
Stricly speaking, the time resolution 6t is finite and the time derivatives in the 
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following text are to be understood as the coarse-grained rates of change of 
the corresponding quantities. However, we shall use the continuum notation 
for simplicity. 

Going over to the interaction representation we get 

Wi{t) = -i[Vi{t),Wi{t)] (15) 
We assume that W{t) factorizes at t = 

W{Q)=p{{))X (16) 

where X is a stationary density matrix describing the internal dynamics of 
the "bath". We again make an idealization: of course, the "bath" is cooling 
as time goes on. However, this cooling is expected to be slow during the 
one-dimensional regime we are interested in and we neglect it altogether ( in 



ref. [13 1 it is found, using the hydrodynamical model, that the temperature 
falls with proper time like r~3 only). Furthermore, the "bath" is not at rest 
with respect to the "system". To deal with this complication we adopt a 
simple ansatz: 

X~ea:p[-^/?P^« + w^)] (17) 

where is the energy- momentum operator of the "bath" and Ui are covariant 
velocities of the outward expansion of the "bath". For ui = U2 = (1, 0) (P!7|) 
becomes the standard canonical density. We assume for definiteness that 
u\ = —u\ in the rest frame of the condensate. Then the exponent becomes 
—P'jHbath ■ The spectrum of the "bath" is red-shifted due to its expansion 
and has an effective temperature T/7. The Lorentz factor 7 is merely a 
phenomenological parameter. 



For A = the problem has the well known solution [T^ 

wP{t) = D[a{t)]\0){0\D%mX (18) 

where 

a{t)=i [' dt'jit')e"''' (19) 
Jo 

is the Fourier component of the classical solution of the equations of motion 
and D{a) is the unitary operator 
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D(a) = exp[aa — aa^] (20) 

Defining 

W{t) = D^[a{t)]Wi{t)D[a{t)] (21) 

and similarly for other operators, we reduce the von Neumann-Liouville equa- 
tion to 

W{t) = -z[V{t),W{t)] (22) 
The perturbative solution of (^2]) is 

W{t) = W{0) + J2{-ir dt,... dtr,[V{h),...,[V{t,,),wm..] (23) 

We work to order O(A^) only and take a trace over the "bath" degrees of 
freedom in order to derive a master equation for the reduced density matrix 
pit) = trf,athW{t) . The perturbative approach is, of course, meaningful only 
when the interaction between the "bath" and the "system" is weak enough 
not to upset the assumed thermal equilibrium of the "bath". The derivation 
starts from the observation that p{t) is obtained from p(0) by the action of 
a linear operator, call it U{t). A simple result is obtained in the so-called 
Markovian approximation , viz. replacing the operator U{t)U^^{t) by its 
value at if: = oo. The main conditions for this to be realistic is that the 
relaxation time of the "bath" is much shorter than that of the "system". 
We adhere to this idealization, although in our problem it is likely to be 
particularly crude, since "bath" oscillations appear slowed down due to its 
expansion. 

Using the notation 

Tnath{...X) = (...) (24) 

we have 

m) = (25) 
since b is linear in 6B. We also assume, for the sake of simplicity, that 
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ibmt')) = 

The master equation is then 



(26) 



p = —i\'^A[d^a, p] + A^K{[a, pa^] + [ap, a^]} + 2Kn\'^[a, [p, a^]] (27) 
where 

POD 

K + iA= / dte'^'{[b{t),b^{0)]) (28) 
Jo 

POD 

Kn = Re dte''^'{b\0)b{t)) (29) 
Jo 

Eq. (^) is formally identical to that governing the behaviour of a damped 
quantum oscillator. Eq. ([19|) implies that generically a(t) becomes indepen- 
dent of time provided t > to , with the rough estimate 

toW ~ 1 (30) 

Thus, for t > to we can drop the tildas in (P7|). In this way we finaly 
reduce our problem to that of a damped oscillator. Of course, we have by no 
means demonstrated that the kinetics of a disoriented chiral condensate is 
equivalent to that of a damped quantum oscillator. We have postulated this 
analogy from the outset, because we feel it may give some insight into the 
real problem. The aim of this section is only to explain the real significance 
of the analogy, by making explicit the dynamical postulates it involves. 

4 THE HIERARCHY OF TIME SCALES 

Let us go to the Schroedinger representation and let us write the reduced 
density matrix in the Glauber representation 

p{t)= Jd'f3\(3)P{(3,M{[3\ (31) 

It is well known that a master equation for p can be converted into a 
Fokker-Planck equation for the Glauber function P. In the case of a damped 
oscillator the solution of the latter equation is known explicitly: 
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P(/?, P,t)=J <fi PiP, tie, I to) (32) 

where 
and 



P(/?,At|e,0 = ^exp[-r/|/3-ee-^(*-*o)p] (33) 

TT 



^-i = ^[l_e-2AM*-*o)] 
= z(cj + A^A) + A^/t (34) 

Assume, that A is so small that at t = to the "system" is still almost coherent, 
viz. P(e,^, to) — (^^(e — tt(to)) . Then the Glauber function at t > to is 
given by the right-hand side of (|33D with ^ = a (to). It is an easy exercise 
to calculate the elements of the density matrix in the occupation number 
representation. One finds 



{N\p\N) = {n + l) 

and 



n 



n+1 



[1 + 0(e-2^ "*)] (35) 



{N\p\M) ~ Ny^M (36) 

We are now in position to identify the different time scales entering our prob- 
lem. First, there is the time to ^ uj~^ required to build up the bremsstrahlung 
field. For a given rapidity cell, the probability of finding quanta is Pois- 
sonian with average N ~ |a(to)P and width Vn. As time goes on the off- 
diagonal terms of the density matrix die out following the law (|36|). The co- 
herence between the relevant distinct multiplicity states of the bremsstrahlung 
field is lost for t > tdecoh, with 

tdecoh ^ [X^k^/n]-' (37) 

Finally, the "system" reaches a stationary state. Indeed, the fluctuation- 
dissipation theorem implies that n^^ = exp{(3uj) — 1. Hence, it follows 
from eqs. (|5D-(|6]) that for 2A^fi:t >> 1 the density matrix is {N\p\M) ~ 
^MN exp{—[3ujN). The thermalization time-scale ttherm is clearly 
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ttherm ^ [A^ft:]"^ (38) 

For large enough multiplicity and small coupling between the "system" and 
the " bath" one finds the following hierarchy of time scales 

U ^ < to < tdecoh << Uherm (39) 

The second inequality is a dynamical constraint. It should be satisfied if the 
transient creation of a coherent classical field is to take place. 

The estimates (|37|) - (|39D have so simple an intuitive meaning that one 
is tempted to believe that their validity transcends the over-simplified model 
used to get them. 



5 SUMMARY AND DISCUSSION 

The foregoing discussion has actually two facets. We have started by try- 
ing to extend the results of to 1+3 dimensions. In this attempt we have 
encountered difficulties, which are at least partly rooted in the assumptions 
we have adopted. Of course, these assumptions may just refiect our prej- 
udices. Anyhow, we have noticed that the difficulties can be circumvented 
if one considers the disoriented chiral condensate as an open quantum sys- 
tem. Hence, we have embarked into the discussion of the time scales entering 
the kinetics of the condensate, insisting on the (presumed) analogy with the 
damped quantum oscillator (or laser under threshold). However, the rel- 
evance of discussing decoherence does not depend on the soundness of our 
original motivation. This question has to be faced in all scenarios. Therefore, 
consider it first. 

As one might expect, the decoherence time scale depends on the (Fourier 
transform of) the linear response function of the "bath" and on the strength 
of the couphng between the "system" and the "bath". A thorough discusion 
of these issues goes beyond the scope of this paper. In the present state of the 
art it would inevitably involve much model building. Let us limit ourselves 
to some rather obvious comments and guesses. 

From the definition of the spectral function k(co') it follows that it satisfies 
the sum rule 
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j dujK{u) = 1 (40) 

Let r denote the size of the bandwidth of the "bath". Roughly speaking 
k{lj) ~ for LV within the bandwidth and zero otherwise. Remember, 
that k{u) has to be calculated setting the temperature of the "bath" to the 
effective value T/^. In a very high-energy collision this effective temperature 
is presumably low and it is perhaps reasonable to guess that the bandwidth 
extends from the lowest frequencies to F, and that the latter is controlled by 
a hadronic scale, say F < 1 GeV . 

In the relativistic regime, the assumed gradient nature of the coupling be- 
tween our "system" and the "bath" imphes that the couphng A = Xqu;, where 
Ao is a dimensionless constant. Since the terms in the action involving the 
field B are supposed to mimic the interactions neglected when one keeps only 
the lowest order term in the chiral action, it is perhaps reasonable to guess 
that Ao ^ J-k/M, where M is the mass scale where the chiral perturbation 
theory breaks down (M ~ 1 GeV, say). Considering a heavy- ion collision 
one should further multiply the coupling term, i.e. the classical current j and 
the coupling A, by the dimensionless constant f^^R. One then finds N (x R^, 
as it should be. 

An extension of the discussion to the non-Abehan case presents additional 
difficulties. However, as long as chiral-symmetry is maintained, the disori- 
entedness of the chiral condensate should not be affected by the decoherence. 

Notice, that at t ~ tdeco/i the multiplicity distribution is still nearly Pois- 
sonian. The coherence between different multiplicity states has been broken, 
but correlations had no time to build up. However, for t > t decoh the evolution 
of the system is no longer governed by the classical equations of motion. 

Let us now return to the scenario outlined in sect. 2. As argued there, the 
classical field should loose its coherence before the end of the one-dimensional 
regime. This would invalidate the argument based on the classical equations 
of motion and maintaining that the pion field will collapse on the sources dur- 
ing the three-dimensional expansion. The condition for that not to happen 
is 

idecoh ~ -R (41) 

This is a stringent condition. It implies that only the quanta belonging to 
the bandwidth of the "bath" may contribute to the observable signature of 
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the disoriented condensate. Thus, the decay products of the condensate are 
expected to be locahzed in rapidity, presumably within an interval a few 
units long. This is to be contrasted with the purely classical result of ref. 
[0], where a rapidity plateau was found 

If our guess for F and Aq is not too far from reality, the condensate 
should decay into a large number of pions, of the order of 10^ for R ^ 3 
fm. Otherwise, the constraint (0) would be hard to meet. An extension of 
the model to i? ^ 1 fm is problematic. If it is nevertheless attempted, then 
for Ao — 0.1 the satisfaction of (|4ll) requires either an unreasonably high 
multiplicity or a very narrow bandwidth. This probably means that one 
should not expect to discover the disoriented chiral condensate in hadron- 
hadron collisions when final states have generic overall topology. Here we 
join the intuition of the authors of refs. P, ||. 

We would like to end this very speculative paper with a few general 
remarks: 

• Contrary to what we have stated at the end of ref. 0, describing 
the creation and evolution of a disoriented chiral condensate in 1+3 
dimensions is not devoid of conceptual challenges. 

• Whatever scenario one is willing to adopt, one should remember that 
the chiral condensate is an open quantum system. The implications of 
the existence of the corresponding hierarchy of time scales are worth 
being examined. 

• The existence of Centauro-like events, if confirmed, would not be a mere 
curiosity. The creation and evolution of a disoriented chiral condensate 
is sensitive to the physics governing the early stage of the collision pro- 
cess and the experimental signature of the phenomenon is particularly 
unambiguous (the neutral/charged ratio). Althouh these events are 
presumably rare, they might convey relatively clean information. 
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